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I. INTRODUCTION: 
WHY FOUR-NEUTRINO MIXING? 

Evidences in favor of neutrino oscillations (see 0|) 
have been found in solar (Homestake Kamiokande 
§, GALLEX [|, SAGE [|, Super- Kamiokande §) and 
atmospheric (Kamiokande, 1MB, Super-Kamiokande, 
Soudan 2, MACROjTj) neutrino experiments and in the 
LSND experiment It is rather well-known that 
these three evidences require the existence of at least 
three small neutrino mass-squared differences (Am^'s): 

Am2^„ - 10-1" - 10-4 cV^ , (1) 
Am2,^ ^10-3- 10-2 eV^, (2) 
Am^g^D'-lO-^-lOeV^. (3) 

This means that at least four light massive neutrinos 
must exist in nature. Following the old principle known 
as Occam's razor ("pluralitas non est ponenda sine ne- 
cessitate" ) we assume the existence of four light massive 
neutrinos, that is the minimal possibility compatible with 
all data. In this case the left-handed components VaL of 
the flavor neutrino fields are superpositions of the left- 
handed components VkL (A; = 1, . . . , 4) of neutrino fields 
with definite mass toj, smaller than a few eV: 

4 

I'aL =^Uakl^kL , (4) 

where U is a 4x4 unitary mixing matrix. 

From the measurement of the invisible decay width of 
the Z-boson it is known that the number of light active 
neutrino flavors is three (see [Tl|]), corresponding to v,,, 
j/^ and (active neutrinos are those taking part to stan- 
dard weak interactions). This implies that the number 
of massive neutrinos is bigger or equal to three. If there 
are four massive neutrinos, in the flavor basis there is one 
sterile neutrino, Vg, that does not take part to standard 
weak interactions. For the flavor index a in Eq. (^) we 
choose the ordering a — e,s, fj,, t. 



Let us emphasize that four-neutrino mixing is very 
interesting both for theory and experiment. From the 
theoretical point of view the existence of a light ster- 
ile neutrino can be explained only with models far be- 
yond the Standard Model and is thus a signal of excit- 
ing new physics. From the experimental point of view 
four-neutrino mixing is attractive because it allows the 
existence of a Am? at the eV scale which generates os- 
cillations in different channels that can be explored with 
high precision in short-baseline experiments. Moreover, 
it allows large CP-violation effects in long-baseline and 
j/-factory experiments (see fl^ ). Four- neutrino mixing 
may also open the possibility for measurable exotic phe- 
nomena, as neutrino decay, neutrino magnetic moment, 
etc. 

From the model-builder point of view, there seems to 
be no difficulty in constructing models with a light sterile 
neutrino (as one can convince oneself with an appropri- 
ate search of the hep-ph electronic archive, resulting in 
too many four-neutrino models to be cited here). How- 
ever, our impression is that these models are constructed 
with assumptions "ad hoc" in order to generate a light 
sterile neutrino, with the possible exception of models 
based on the existence of a mirror world (see jl^]) and 
models based on the existence of large extra dimensions 
(see [0), in which light neutral fermions that mix with 
the ordinary neutrinos seem to be naturally allowed. 

Our approach is phenomenological: we extract infor- 
mation on four-neutrino mixing (masses and mixing an- 
gles) from the available data. 

II. ALLOWED FOUR-NEUTRINO SCHEMES 

The six types of four-neutrino mass spectra with three 
different scales of Am? that can accommodate the hier- 
archy Am2^„ < Aml,^^ « Amls^^, (see Eqs. 
are shown qualitatively in Fig. nl In all these mass spec- 
tra there are two groups of close masses separated by 
the "LSND gap" of the order of 1 eV. The six schemes 
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are divided into four schemes of class 1 (I-IV) in which 
there is a group of three masses separated from an iso- 
lated mass by the LSND gap, and two schemes of class 
2 (A, B) in which there are two couples of close masses 
separated by the LSND gap. In each scheme the smallest 



mass-squared difference (Am^^- = — m|) corresponds 
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served in the LSND experiment ( Amli 
in the schemes of class 1, Ami 
Ato|2 in the schemes of class 2). 

It has been show that the four schemes of class 1 are 
disfavored by the data if also the negative results of short- 
baseline accelerator and reactor disappearance neutrino 
oscillation experiments are taken into account p5| p7|,p[ . 
The reason is simple: in order to explain solar and at- 
mospheric data with neutrino oscillations Ve and Vfj, must 
have large mixing with the group of three neutrinos and, 
consequently, small mixing with the isolated neutrino. 
In this case, short-baseline oscillations occurring through 
the mass-squared difference Ato^j^ , i. e. through the inter- 
ference of the wave functions of the isolated neutrino and 
the three grouped neutrinos, are strongly suppressed. It 
turns out that the negative results of short-baseline accel- 
erator and reactor disappearance experiments imply that 
the amplitude of short-baseline i^^ transitions is 

suppressed below the value measured by LSND ||l^,|l^,^ . 
Hence, the schemes of class 1 cannot explain all neutrino 
oscillation data. 

On the other hand, the two four-ncutrino schemes of 
class 2 are compatible with the results of all neutrino os- 
cillation experiments if the mixing of with Pi and V2, 
the two mass eigenstates responsible for the oscillations 
of solar neutrinos, and the mixing of i'^ with and 1^4, 
the two mass eigenstates responsible for the oscillations 
of atmospheric neutrinos, are large Hence, we 

have the two possibilities illustrated in Fig. ^, in which 
Ue and Ufj, are depicted close to the two massive neutri- 
nos with which they have large mixing and it is explicitly 
indicated that the mixing of Vr and Vs is unknown. In- 
formation on the mixing of Vt and Vs, as well as more 
detailed information on the mixing of and v^, can be 



trino data, and, in the future, of the data of long-baseline 
experiments. 

The appropriate formalism for the treatment of solar 
and atmospheric neutrino oscillations, including matter 
effects, in the framework of the two four-neutrino schemes 
A and B has been presented in Ref. ||l8|. In the next 
section we present the results of a fit of solar neutrino 
data jl9). 



III. FIT OF SOLAR NEUTRINO DATA 

In the framework of the two four-ncutrino schemes A 
and B, solar neutrino oscillations depend on only one 
mass-squared difference, Ato^^j^, but there are several 
degrees of freedom coming from the 4x4 mixing ma- 
trix, that can be parameterized in terms of six mixing 
angles and three CP-violating phases (three additional 
CP-violating phases are possible if massive neutrinos are 
Majorana particles, but can be neglected in the present 
contest because they do not have any effect in neutrino 
oscillations). However, the data of solar neutrino experi- 
ments together with those of short-baseline disappear- 
ance experiments indicate that the elements Ue3 and Uei 
of the neutrino mixing matrix are very small (l^j2^ . This 
implies that two mixing angles are small and can be ne- 
glected in the study of solar (and atmospheric) neutrino 
oscillations [|l8|. Furthermore, the CP-violating phases 
can be neglected because their effects are washed out 
by the average over neutrino energy and source-detector 
distance. Therefore, we consider the mixing matrix U in 
Fig. I ||l8|,|9|, where 1^12, 1^23, 1^24, '^34 are four mixing 
angles and Cy = cos z9y and = sin t9.y . The rows and 
columns of the mixing matrix U correspond, respectively, 
to the neutrino flavors I'e , i^s , t'^j , i^t and to the mass eigen- 
states 1^1 , 1^2 7 f 3 , 1^4 • 

Since solar neutrino oscillations are generated by the 
mass-square difference between V2 and 1^1 and Uei = 
cost9i2, Ue2 = sin'(?i2, the survival of solar Ve^s, P^^-,^^, 
mainly depends on the mixing angle z?i2, with the only 
correction due to four-neutrino mixing in the matter po- 
tential that is given by 



V = Vcc + COS^ l923 COS^ §24 Vl 



NC ; 



(5) 



obtained through the fit of solar and atmospheric neu- where Vcc — V^GpNg and Vnc = — a/2G_fA^„/2 are 
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Experiment 


Rate 


Units 


BP98 


Homestake 


2.56 ± 0.23 1 




SNU 


7.8 ±1.1 


GALLEX+SAGE 


72.3 ±5.6 g 




SNU 


130 ±7 


Kamiokande 


2.80 ±0.38 [ 




10*^ cm-^s"^ 


5.2 ±0.9 


S up er- Kamiokande 


2.45 ±0.08 [ 




10"^ cm-^s-i 


5.2 ±0.9 
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the usual charged-current and neutral-current potentials 
(see §). 

The mixing angles i?23 and 1^24 determine the relative 
amount of transitions into sterile Vg or active and I'r 
(that cannot be distinguished in solar neutrino exper- 
iments, because their matter potential and their inter- 
action in the detectors are equal, due only to neutral- 
current weak interactions). The active/sterile ratio and 
solar neutrino oscillations in general do not depend on 
the mixing angle 'i?34, that contribute only to the dif- 
ferent mixings of and Vt. Indeed, the mixing of I's 
with vi and i>2 depends only on t?i2 and the product 
cos 7?23 cos 1^24. Moreover, instead of and Vt, one can 
consider the linear combinations 

[ ^<^ \ ^ [ -sini9 -cost) \ f S34 C34 W ^ 
\vb ) \ cos.d -sm{} ) \ C34 -S34 ) \i^T ) ' 

with tani9 = sin 1^24/ tan 1^23 • The mixing of Va and i/fc 
with vi and V2 is given by Uai = -S12 ^1 ^ ^23034, Ua2 = 
ci2-\/l — C23C24, Ubi = Ub2 = 0. Therefore, the oscilla- 
tions of solar neutrinos depend only on i?i2 and the prod- 
uct COSI?23 COSl?24. If < COS ??23 COS Z?24 < 1, Solar Z^e's 

can transform simultaneously into sterile neutrinos and 
in the linear combination i/^ of active z/^ and Vr , with the 
limiting cases of pure two-generation i/^ — > Va transitions 
if cos z923 cos i?24 = {Usi = Us2 = 0, Uai = -sinz9i2, 
Ua2 — cos 1912), and pure two-generation Ve Vs transi- 
tions if cos i?23 C0S1?24 = 1 {Usl = — sin 7^12, Us2 = C0S??12 

and Uai = Ua2 = 0). 

Since the mixing of v^. with vi and V2 is equal to the 
one in the case of two-generations (with the mixing angle 
i9i2), the mixing of with vi and V2 is equal to the one 
in the case of two-generations times cos i?23 cos i?24 and 
the mixing of Va with vi and V2 is equal to the one in the 
case of two-generations times \/l — cos^ 1^23 cos^ i?24, it is 
clear that in the general case of simultaneous Ve — >■ 
and Ve — > Va oscillations the corresponding transition 
probabilities are given by 

P,^^,^ = cos" l923 cos" ^24 (1 - Pu^-.u^ , (6) 
P^^_^^ = (1 - cos" 7923 cos" d2A) (1 " P^.-^^J • (7) 

These expressions satisfy the relation of probability con- 
servation P^^^i,^ + Pi/^^jy^ + Pv^-ty^ — 1. We cal- 
culated the survival probability P^^^v^ in the range 
10~"eV" < Am" < 10"" eV", 10""* < tan" 1^12 < 10 
and < cos" ^23 cos" '!924 < 1, taking into account matter 



effects for 10"^ eV" < Ato^^ < 10"^ eV" and the possi- 
ble regeneration of I'e's when the flux of solar neutrinos 
crosses the Earth (for details see Ref. jigf). 

In order to calculate the allowed regions for the param- 
eters Am2i, tan" '{)i2 and CI3CI4, we have used data on the 
total event rates measured in the Chlorine experiment at 
Homestake [|| , in the two Gallium experiments GALLEX 
[Q and SAGE and in the water Cherenkov detectors 
Kamiokande and Super-Kamiokande Q shown in Ta- 
ble |. We used also the zenith angle distribution of the 
events and the electron recoil energy spectrum measured 
in the Super-Kamiokande experiment (825-day data sam- 
ple) Ql . For the calculation of the theoretical expecta- 
tions we use the BP98 standard solar model of Ref. . 
The calculation of the theoretical rates, the statistical 
treatment of the data and the results of different fits of 
rates only, rates and zenith angle distribution, rates and 
recoil electron energy spectrum are described in Ref. [ p^ . 
Here we present only the result of the global x" fit of all 
data. 

Figure || shows the allowed regions at 90% and 
99% CL in the tan" i?i2-Am2i plane for C23C24 — 
0,0.2,0.4,0.5,0.71, that we calculated in the following 
way. We computed for a grid of points in the three- 
dimensional tan" i?i2~Am|2"c"3c"4 parameter space the 
expected values of the observables and with those and 
the corresponding uncertainties we constructed the func- 
tion x"(tan" i?i2, Am"2, C23CI4). We found its minimum, 
Xmin, in the full three-dimensional space considering as 
a unique framework both MSW and vacuum oscillations. 
The allowed three-dimensional regions for a given CL are 
then defined as the set of points satisfying the condition 

x"(tan"i?i2,Am?2,c^3C^4)-xLn < Ax'(CL,3dof) (8) 

where Ax"(CL,3dof) = 6.25,11.36 for CL=90%, 99%, 
respectively. In Fig. ^ we plotted the sections of such 
volumes in the tan" z?i2-Ato2i plane for six selected val- 
ues of C23C24. 

The global minimum of x", Xmin — 28.8 with 23 
degrees of freedom corresponding to 28 data points (4 
rates plus 6 Super-Kamiokande zenith-angle bins plus 
18 Super-Kamiokande energy spectrum bins) minus 2 
normalization factors (for the Super-Kamiokande zenith- 
angle and energy spectrum data) minus 3 fitted param- 
eters, lies in the LMA region for cos" 1^23 cos" 1^24 = 0, 
corresponding to pure two-generation v,, — > Ua oscilla- 
tions, and tan"i9i2 = 0.35, Am^^ = 2.3 x 10"^ eV". 
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However, as one can see from Fig ||, the of the 
best-fit point in the LMA, LOW and VO regions rises 
steeply when cos^ i?23 cos^ 1^24 is increased. Therefore, 
the 99% CL LMA, LOW and VO regions disappear 



for 



-23 



close to one, 



0.78, 0.70, 0.83, respectively. For values of 
corresponding to almost pure two- 
generations Ve — > i^s transitions, only the SMA region is 
allowed. This is expected, because only the SMA solu- 
tion is allowed in the two-generation i^s fit of solar 
data (see [^). Notice, however, that the statistical anal- 
ysis is different: in the two-generation picture the pure 
j^e— ^active and I'e^sterile cases are analyzed separately, 
whereas in the four-neutrino picture they are taken into 
account simultaneously in a consistent scheme that al- 
lows to calculate the allowed regions with the prescription 
given in Eq. (||). We think that the agreement between 
the results of the analyses with two and four neutrinos 
indicate that the physical conclusions are quite robust. 

The disappearance of the LMA and LOW regions for 
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close to one is due to the following reason. 



Unlike active neutrinos which lead to events in the water 
Cherenkov detectors by interacting via neutral current 
with the electrons, sterile neutrinos do not contribute 
to the Kamiokande and Super-Kamiokande event rates. 
Therefore a larger survival probability for neutrinos 
is needed to accommodate the measured rates. As a con- 
sequence a larger contribution from ®B neutrinos to the 
Chlorine and Gallium experiments is expected, so that 
the small measured rate in Chlorine can only be accom- 
modated if less ^Be neutrinos are present in the flux. 
This is only possible in the SMA solution region, since in 
the LMA and LOW regions the suppression of ''Be neu- 
trinos is not enough (the survival probability of VeS as a 
function of neutrino energy is rather flat). 

Figure | shows that the LMA region at 99% CL ex- 
tends to high values of Am^i, even above 10~^eV^ for 
CI3C24 ^ 0.1. Since the atmospheric mass squared differ- 
ence Ato^jj^ hes between 10~^ and 10~^ eV^ (see 0), one 
may wonder if the solar and atmospheric mass squared 
differences may coincide and three massive neutrinos may 
be enough for the explanation of solar, atmospheric and 
LSND data. The answer to this question is negative, be- 
cause in the high-Amli part of the 99% CL LMA region 
the mixing angle 6*21 is large, 0.3 < sin^(02i) ^ 0.7, and 
in this case disappearance of v^s should be observed in 
long-baseline reactor experiments, contrary to results of 
the CHOOZ experiment. In other words, the results 
of the CHOOZ experiment, that have not been taken into 
account in the present analysis, forbid the part of the 99% 
CL LMA region that extends above Amji ~ 10^"^ eV^. 
For this reason we cut the plots at this value. 
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IV. CONCLUSIONS 

We have considered the two four-neutrino schemes A 
and B in Fig. ^ that are compatible with the results of 
all neutrino oscillation experiments. These two schemes 
are completely equivalent for neutrino oscillation exper- 
iments, but have different phenomenology in /3 and neu- 
trinoless /3/3 decay experiments (see p], p4^ , p5| ). 

In general the four-neutrino schemes A and B allow si- 
multaneous transitions of solar z/e's into active I'^'s, z^r's 
and sterile Vs, as well as simultaneous transitions of at- 
mospheric I'p's into active Vr^s and sterile i^s's. 

We have fitted the data of solar neutrino experiments 
in terms of neutrino oscillations, that in the schemes un- 
der consideration depend on three parameters, Am^i, 
tan^ z9i2 and cos^ 'i?23 cos^ t?24- The allowed regions of 
these three parameters are depicted in Fig. |4|, from which 
one can see that the SMA region is valid for all values of 
cos^ 'i?23 cos^ 'i?24 («-e. any combination of i/e^active and 
z^e^sterile transitions), whereas the LMA, LOW and VO 
solutions disappear for values of cos^ cos^ close to 
one, where Ve — *■ f^s transitions are dominant. 

Other authors have recently performed fits of the atmo- 
spheric neutrino data in the framework of the two four- 
neutrino schemes A and B [^,^ . We expect that in the 
future a combined fit of solar and atmospheric neutrino 
data will allow to constraint further the mixing of four 
neutrinos. 
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